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Abstract 

We study the phase diagram of Ri_ x A x Mn03 (R=La, Pr, Nd, Sm ; A=Ca, Sr, 
Ba) by taking into account the degeneracy of e g orbitals and the anisotropy of 
the transfer integral. The electron-electron interaction is treated in the mean 
field approximation with the optimization of the spin and orbital structures. 
The global phase diagram is understood in terms of the two interactions, i.e., 
the super exchange interaction for small x and the double exchange interac- 
tion for larger x modified by the orbital degeneracy. The dimensionality of 
the electronic energy band resulting from the orbital structure is essential to 
determine the phase diagram. The effects of the Jahn- Teller distortion are 
also studied. 
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The understanding of the rich phase diagram in Ri-^A^MnOa (R=La, Pr, Nd, Sm ; 
A=Ca, Sr, Ba) is indispensable for the discussion of its physical properties including the 
colossal magnetoresistance. These materials have been considered to be the model system 
of the double exchange mechanism 0-0, i.e., the t2 9 spins are aligned parallel in order to 
minimize the kinetic energy of e g electrons, which are strongly Hund coupled to t 2g spins. 
However this simple picture of the ferromagnetic phase has been questioned recently by 
several authors 0-^4], who stress the importance of the other interactions in addition to 



the double exchange one. The parent compound LaMnOs is an insulator with the A-type 
antiferromagnetic ( AF) ordering and the Jahn- Teller ( JT) distortion [|l^,|n| , while it should 
be metallic when only the Hund coupling is considered. Roughly speaking there are two 
streams of thinking on this issue. One is to regard the JT distortion to be of the primary 
importance MM, which removes the degeneracy of the orbitals. The other is to stress the 



strong correlation effects of the e g electrons [11 14]. In this picture the parent material is 



regarded as a Mott insulator, and the effective Hamiltonian is derived to study the spin 
and orbital structures JT^JT^] • In the local density approximation (LDA) and LDA+U band 



calculations for a; = ||, where the effect of the electron correlation is included in a kind of 
mean field approximation, it is concluded that the JT distortion of the (3x 2 — r 2 ) / (3y 2 — r 2 )- 
type is important for the A-type spin structure observed experimentally. On the other 
hand, a recent exact diagonalization study of the effective Hamiltonian [|17| suggests that 
the correlation of the (3x 2 — r 2 ) / (3y 2 — r 2 )-type or (z 2 — x 2 ) / (y 2 — z 2 )-type orbital structure 
becomes remarkable in the A-type AF phase. As for the doped case (x ^ 0), the system 
becomes ferromagnetic metal for x > 0.175 The simple double exchange mechanism 

is considerably modified as shown below when the anisotropy of the transfer integrals 
between the e g orbitals is taken into account. Especially it has been suggested that the orbital 



degrees of freedom might remain disordered down to low temperatures [0 to explain the 
anomalous physical properties. Hence the origin of the ferromagnetism should be reexamined 
taking into account the orbital degeneracy. Near x = 0.5 the charge ordering accompanied 
with the spin and orbital orderings has been observed. With x increased further (x = 0.6), 



the A-type AF structure again appears, which shows quasi-two dimensional metallic behavior 



In this paper we present an extensive study on the phase diagram of Ri-^A^MnOs 
(R=La, Pr, Nd, Sm ; A=Ca, Sr, Ba) in the mean field approximation, which treats both the 
super and double exchange interactions in a unified fashion at zero temperature. Because 
the Coulomb interaction is the largest interaction and also the JT distortion disappears for 



x > 0.15 pOfl , we first study the model with only electron-electron interactions. The spin 
and orbital structures are optimized, and the global phase diagram is given in the plane of 
x ( the concentration of the holes ) and J s ( the super exchange interaction between the 
t 2g spins ). The effect of the JT distortion, which turns to be important for x = 0, is also 
studied. 

We set up the three-dimensional cubic lattice consisting of the manganese ions. Two 
kinds of the e g orbital( 7,7') are introduced on each site, and the t 2g electrons are treated 



as a localized spin with S = 3/2. The Hamiltonian without the JT coupling is given by |T2 

H = E (tl/ 'AvAm* + h - c ) 

<ij>,a,-y,7' 

+ U E ^ 7 T n n! + U ' E n ia. n ib + 1 E d la*4bv> d iaa> d ib* 

+ Jh Y,s 1 -s^ + j s E s^-s?° . (i) 

i <ij> 

d\ la is the operator which creates an electron with spin cr(=T, I) in orbital 7(= a, b) at 
site i, and Si is the spin operator for the e g electron defined by Si = \ *}2ao''yd\ (T1 (j (Ta idi CT i 1 . 
The electron transfer integral t]J , which is estimated by considering the oxygen 2p orbitals 
between the nearest Mn-Mn pair, is represented as c]J to, where c[J is the numerical factor 
depending on the orbitals and to is estimated to be 0.72eV which we choose the unit of energy 
below (to = 1) ||12|| . The second line shows the electron-electron interaction terms where 
U, U' and / is the intra-, inter-orbital Coulomb interactions, and inter-orbital exchange 
interaction, respectively. This interaction can be rewritten as -aEi(Si(r) + 2^^ 2s ( T )) — 
P Hi Ti{r) 2 |L4j] . Here the spin operator Si and the iso-spin operator Tj = \ D 77 'o- dl^a^ydiay 



3 



for the orbital degrees of freedom are introduced, and the two positive coefficients a and j3, 
which are defined by a = 2U/3 + U'/3 — 1/6 and (3 = U' — 1/2 , represent the interaction 
to induce the spin and iso-spin moments, respectively. The last line is the sum of the 
Hund coupling and the AF interaction between the nearest neighboring t2 g spins. Here we 
adopt the mean field approximation by introducing the order parameters (Si), (S^ 29 ), and 
(Tj). These order parameters are determined to optimize the mean field energy at zero 
temperature. For both spin and orbital, the four types of the ordering are considered, that 
is, the ferromagnetic (F-type) ordering, where the order parameters are uniform, and the 
three AF-like orderings, i.e., the layer- type (A-type), the rod-type (C-type) and the NaCl- 
type (G-type) AF orderings. Hereafter, types of the orderings are termed as, for example, 
(spin:C), and so on. 

In Fig. 1, the spin and orbital phase diagram is shown for the set of parameters a = 
70 >> (3 = 2.5. In this rather extreme case the regions dominated by the super exchange 
and double exchange interactions are separated, and it is easy to obtain the physical picture. 
a/ (3 » 1 corresponds to the situation where 1/{U' — I) » 1/U, l/(U' + 1). In Fig. 1, the 
spin structure changes as F — > A — > C — > G, as J s increases, which is consistent with the exact 
diagonalization study [T7|]. We begin with the discussion of the parent material (x = 0) where 



only the super exchange interaction is effective. For spin:A, which is observed in LaMnOs, 
the most stable orbital structure is the orbitahG {[3z 2 — r 2 ] + [x 2 — y 2 ])/([3z 2 — r 2 ] — [x 2 — y 2 }) 
as shown in Fig. 1. When the ratio a/ (3 is changed, the orbital changes continuously 
from orbitahG ([3z 2 - r 2 } + [x 2 - y 2 })/([3z 2 - r 2 } - [x 2 - y 2 }) for a/(3 » 1 to orbitahG 
(y 2 — z 2 )/(z 2 — x 2 ) for a/ (3 ~ 1, and to {3z 2 — r 2 ) for a/ [3 = 0. For the actual compound, 
we expect a>j (3 and orbital:G(?/ 2 — z 2 )/(z 2 — x 2 ) is the most stable. The experimentally 
observed orbital:G(3a; 2 — r 2 )/(3y 2 — r 2 ) is never the most stable solution, which can be 
understood as follows. There are three possibilities for the intermediate states of the super 
exchange processes, i.e., the occupancy of the two orbitals (a) with the parallel spins (the 
energy U' — I) or (b) antiparallel spins (U' + /), and (c) the double occupancies of the same 
orbital (U) jl2|,[l3l,pl| , ^2l . Then the relative importance of the states (a) compared with (b) 



and (c) is changed when a/ (3 is changed. Let us compare the energy gains due to the super 
exchange processes in orbital:G (y 2 — z 2 ) / (z 2 — x 2 ) and orbital:G (3x 2 — r 2 ) / (3y 2 — r 2 ). For 
the processes using states (a) and (b), the magnitudes of the transfer integrals and hence 
the energy gain are the same, while for the process using (c), the energy gain is always larger 
for (y 2 — z 2 )/(z 2 — x 2 ) compared with (3x 2 — r 2 )/(3y 2 — r 2 ). Then there is no chance for 
{3x 2 — r 2 )/{3y 2 — r 2 ) to be the most stable structure for any value of a/ (3. Hence the JT 
coupling is important in addition to the electron-electron interactions at x = 0. We introduce 
the JT distortion observed experimentally and its coupling to the e g electrons. We consider 
the two shorter Mn-0 bonds and the one longer bond in the MnC>6 octahedron, and its bond 
length is represented as di ong = d (l + 0.056) and d s h ort = d (l — 0.028), respectively, as we 
follow the structural data [15] . The change of the transfer integrals is estimated in terms of 



the dependence of the 3d — 2p hopping t p d on the distance d as t p d oc <i~ 7 / 2 . We vary 
the splitting of the energies (g) between the two orbitals as the parameter and found 
[1] The wave functions are saturated to become (3x 2 — r 2 ) / (3y 2 — r 2 ) when g is about the half 
of the transfer energy t - This value is much smaller than what is expected in the absence 
of the electron-electron interactions. The magnitude of the isospin moment |Tj| is already 
induced almost fully in terms of the strong electron-electron interactions, and the role of the 
JT coupling is to fix the direction of Tj. [2] The spin: A is stabilized relative to spin:F by JT 
distortion. The phase boundary J S (FA) between A and F is shifted from J S (FA) = 0.014 for 
g = to J s = 0.007 for g = 1.0. This tendency is in agreement with the other calculations 



PJT0||, but the physics is different. In the band calculation, the ground state without the 
JT distortion is the ferromagnetic metal and the enhanced AF exchange between layers is 
due to the reduced double exchange interaction by JT distortion 0. In our calculation, on 
the other hand, only the super exchange interactions are relevant because the large gap has 
been already opened up due to the strong electron-electron interactions. The stabilization 
of the A-AF is clearly understood in terms of the change of the super exchange interaction 
by the orbital rearrangement. 

Now let us turn to the doped case (x ^ 0). The orbital structure in spin:F is quite 



sensitive to the carrier concentration, that is, it changes continuously as x increases from 
orbital:G(x 2 — y 2 )/(3z 2 — r 2 ) near x = to (x 2 — y 2 ) for x = 0.3, and to orbital:A([3,2 2 — 
r 2 ] + [x 2 - y 2 ])/([3z 2 - r 2 } - [a; 2 - y 2 }) for 0.3 < x < 0.8 and finally (3^ 2 - r 2 ) for x = 0.8, 
as shown in Fig.l. On the other hand the orbital in spin: A and spin:C almost remains 
(x 2 — y 2 ) and (3,2 2 — r 2 ), respectively, expect for x = 0, in contrast to spin:F case. The 
phase boundary J S (FA) increases linearly near x = 0, and turns to decrease to have a 
minimum J S (FA) = at around x = 0.3, where both spin:F and A have the (x 2 — y 2 ) 
orbital. The linear increase is due to the difference in the location of the band edges for 
spin:F and spin:A structures. This feature remains true even when the canting in the 
spin:A is taken into account because it gives the energy gain only of the order of ~ x 2 
||. The minimum of J S (FA) around x = 0.3, separates rather clearly the two regions 
dominated by the super exchange (x 0.3) and the double exchange interactions (x ^ 0.3). 
In the doped case the shape of the density of states and the Fermi energy is crucial to 
determine the double exchange energy, which depends on both spin and orbital structures. 
Especially The dimensionality and the van- Hove singularities of the density of states depends 
strongly on the orbital structure. Therefore, in that sense, the double exchange mechanism 
is considerably modified from the conventional one, when the anisotropy of the transfer 
integrals and the electron-electron interaction are taken into account. In the region of 
x < 0.3, the orbital:G(x 2 — y 2 )/(3z 2 — r 2 )- and orbital:F(x 2 — y 2 )-type structures realized 
in the spin:F and spin:A phases respectively, and the band becomes two-dimensional. Here, 
the density of states has a logarithmic singularity at e = and decreases monotonously 
as \e\ with the steps at e — ±3to- In the low carrier concentration region, the above type 
of the density of states is more favorable than the three dimensional one due to the step 
at the band edge. In this case there is no difference in the kinetic energy of the carriers 
between spin:F and spin:A, and J s favors spin:A. Then the F region in Fig. 1 for x < 0.3 is 
stabilized by the super exchange interaction. In the region of x > 0.3, the orbital structure 
is rearranged as orbital:A([3^ 2 — r 2 ] + [x 2 — y 2 ])/([3z 2 — r 2 } — [x 2 — y 2 ])-type, where the 
band structure is essentially three dimensional, but the density of states has two peaks at 

6 



e = ±2to and resembles that of the one- dimensional band. Eventually the (3z 2 — r 2 ) orbital 
appears at x — 0.8 and gives the one dimensional- like band along the z-axis where the 
density of states has the peak structures at e — ±3to- Then the adjusting of the orbital 
structure with increasing x occurs in order to minimize the kinetic energy, i.e., the center of 
mass for the occupied states. When one consider the occupied orbital, the energy position 
of the band edge does not depend on the dimensionality of the dispersion. Then only the 
shape of the density of states matters, and one dimensional-like dispersion is advantageous 
in this viewpoint. For spin:C, the orbital is (3z 2 — r 2 ) almost always except at very small 
x. This can be easily understood because the spin structure allows the electron motion only 
along the z-axis. For spin:G, the electron motion is blocked in all directions and the double 
exchange energy gain is absent. Then the electronic energy does not depend on the orbital 
structure in the limit of strong electron-electron interaction. 

In Fig. 2, we present the calculated phase diagram in the case of a = 8.1, [3 = 2.5, which 
is more relevant to the actual manganese oxides. This set of parameters is complementary 
to that in Fig. 1 because both a/ (3 and a are smaller. In comparison with the results in 
Fig. 1, the spin:F region dominated by the super exchange interaction, is extended to the 
region with larger x. This results in the merging of the super exchange and double exchange 
regions. First consider the spin:F state. At x = 0.0, the orbital structure in the spin:F 
phase is the same as that in Fig. 1. What is different from Fig. 1 is that as x is increased 
the orbital becomes orbital:C i\3z 2 — r 2 ] + [x 2 — y 2 ])/([3z 2 — r 2 } — [x 2 — y 2 ])-like structure 
instead of (x 2 — y 2 ) and to the orbital: A ([3z 2 — r 2 } + [x 2 — y 2 ])/([3z 2 — r 2 ] — [x 2 — y 2 ]) around 
x = 0.6. The two dimensional (x 2 — y 2 )-like structure appears for x = 0.9 where the mixing 
of the two bands is important. 

As for the spin:A, C, G, the orbital remains basically the same as in Fig. 1. Then it is 
concluded that the orbital structure is sensitive to the interactions only in spin:F. This is 
related to the the degeneracy of the orbital structures ]14|]. Actually all the obtained orbital 
structures in Figs. 1 and 2 belong to the lowest degenerate states discussed in [III . 



We now discuss the comparison between the mean field phase diagram in Fig. 2 and 



the experiments. In the present mean field calculation, the ferromagnetic phase, i.e. the 
spin:F phase is growing up with increasing x from the insulating state, and it becomes most 
remarkable around x = 0.3, as shown in Fig. 2. The global feature of the spin:F phase is 
quite consistent with the experimental results in Lai_ x Sr x Mn03, where the ferromagnetic 
phase appears at about x = 0.08 and it survives up to x = 0.5. It is worth to note that, 
however, the origin of the ferromagnetic phase is far from the conventional double exchange 
mechanism, i.e., both the super exchange and the double exchange interactions modified 
by the orbitals are relevant in the region 0.2 < x < 0.4. The orbital ordering in spin:F, if 
observed experimentally by neutron and/or X-ray diffraction, will give important clues to 
the interactions because it depends sensitively on the parameters as described above. The 
another possibility is the orbital fluctuation is so large that the orbital liquid state is realized 



T4j. The RPA analysis of the mean field solutions is left for the future work. 

The another implication to the experimental results is about spin:A phase appearing 
around x > 0.5. In Ndi_ x Sr x Mn03, the ferromagnetic metallic phase is realized up to about 
x = 0.48 and the CE-type AF structure with the charge ordering tunes up [H|r9]. With 



further increasing of x, the metallic state with spin:A again appears at about x=0.53, and 
the large anisotropy in the electrical resistivity is observed in this phase. The similar metallic 
phase accompanied with spin: A is also reported in Pri_ x Sr x Mn03 |19|]. Although the charge 



ordered phase is not considered, i.e., the long range Coulomb interaction is neglected, in 
the present work, the global phase change, as spin:A insulator (x ~ 0) — > spin:F metal 
(0.1 i^xi^O.5) — > spin:A metal (0.55 <,x), is well reproduced when we fix J s to be around 
0.02eV, which is a reasonable value, as shown by the broken line in Fig. 2. It is predicted 
that the (x 2 — y 2 )-type orbital structure should be realized in this spin:A metallic phase, in 
contrast to the insulating phase with spin:A appearing at x — 0.0. Furthermore, the AF 
interaction between layers is expected to be enhanced in comparison with that at x = 0.0, 
because the ferromagnetic interactions originated from the both double exchange and super 
exchange interactions are prohibited in this direction. By the same reason it is predicted 
that the spin canting does not occur in this spin: A metal as observed in fli~9 |, in contrast to 



S 



the spin:A in the small x region which has been discussed by de Gennes 

In summary we have studied the phase diagram of Lai-^Sr^MnOs in the plane of x ( hole 
concentration ) and J s ( AF exchange interaction between the t^g spins) in the mean field 
approximation. The global features can be understood in terms of the interplay between the 
super exchange and the double exchange interactions which are considerably modified with 
taking the orbital degrees of freedom into account. The dimensionality of the energy band 
attributed to the orbital structure plays essential roles to determine the phase diagram. The 
orbital structure is sensitive to changes of the carrier concentration and the interactions only 
in the ferromagnetic state, which suggests the importance of the ordering/disordering of the 
orbital on the origin of the ferromagnetism in the perovskite manganites. 
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Figure captions 

Figure 1. The mean field phase diagram in the plane of the carrier concentration (x) and 
the antiferromagnetic interaction J s between the ti g spins. The strength of the interactions 
are set as a = 70 >> (3 = 2.5. The schematic orbital structure in the each phase is also 
shown. 

Figure 2. The calculated mean field phase diagram with a = 8.1 and j3 = 2.5 case. 
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